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Abstract 

A spin-fermion model relevant for the description of cuprates ladders is studied in a path integral 
formalism, where, after integrating out the fermions, an effective action for the spins in term of 
a Fermi-determinant results. The determinant can be evaluated in the long-wavelength, low- 
frequency limit to all orders in the coupling constant, leading to a non- linear a model with doping 
dependent coupling constants. An explicit evaluation shows, that the spin-gap diminishes upon 
doping as opposed to previous mean-field treatments. 

PACS numbers: 71.10.Fd,71.27+a,75.10.Jm 
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I. INTRODUCTION 



Doped quantum antiferromagnets (QAFM) constitute a major unresolved problem in 
condensed matter physics, which is at the center of current research since the discovery 
of high T c superconductivity! . In particular, the case of a doped spin liquid -where no 
symmetry is spontaneously broken- is very challenging, since the starting point, the spin 
liquid state, cannot be described by a classical Neel state. 

This problem is not only of theoretical relevance. CU2O3 ladders are present in 
Sri4_ a; Ca a ;Cu2404i and many experiments support the presence of a spin gap and a finite 
correlation lengthiHilii'B, two crucial ingredients signaling a spin liquid state. With isova- 
lent Ca 2+ substitution of Sr 2+ holes are transferred from the CUO3 chains to the ladders!, 
increasing the conductivity of the latter. The spin gap, as measured by Knight shift or 
NMR experimentsi'&i is seen to diminish. With increasing doping, superconductivity is ul- 
timately stabilized under pressureill^, a phenomenon that suffices to justify the interest for 
the subject. 

The simplest model which is believed to grasp the physics of the problem is the t — J 
model on a two leg ladder. It is believed in general that this system evolves continuously 
from the isotropic case to the limit of strong rung interaction. In this limit some simplifying 
pictures are at hand: without doping the gap is the energy of promoting a singlet rung to 
a triplet (~ J±). Interaction among the rungs leads eventually to the usual magnon band. 
Upon doping the systems shows two different kinds of spin excitations00. One is still the 
singlet-triplet transition as before, the other kind corresponds to the splitting of a hole pair 
into a couple of quasiparticles (formed by a spinon and an holon), each carrying charge 
+ |e| and spin 1/2. The number of possible excitations is proportional to (1 — 5) (for the 
magnons) and 5 (for the quasiparticles), respectively, where 8 is the number of holes per 
copper sites. For this reason, at low doping concentration, the magnon gap will be the most 
important in influencing the form of the static susceptibility or dynamical structure factor. 

First Sigrist et. al. and more recently Lee et. al. attacked the problem ultimately 
with some sort of mean field decoupling. Their results agree in predicting an increase of the 
magnon gap (Am, originated from the singlet-triplet transition), while Lee et. al. were also 
able to calculate a decrease of the quasiparticle gap (Aq P originated from the splitting of a 
hole pair) for small doping concentrations. 
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In contrast to the mean-field results above, Amnion et. al. obtained a decrease of the 
magnon gap and an almost doping independent Aqp using temperature density matrix 
renormalization group (TDMRG). As already mentioned, a decrease of the spin gap is also 
observed in a number of experimentsi'SH. 

In this paper we concentrate on the behavior of the magnon gap upon doping. Due to 
the contradiction above it is imperative to go beyond mean field and include the role of 
fluctuations in a controlled manner. A mapping from an AFM Heisenberg model to an 
effective field theory, the non linear a model (NLcrM), proved very efficient in describing 
the magnetic properties of two dimensional spin lattices^, chains^, and ladders^. This 
mapping was extended in Ref.0 to the case of a doped two dimensional QAFM using a 
procedure that we will closely follow. 



II. MAPPING TO AN EFFECTIVE SPIN ACTION 

Since no satisfactory analytical treatment of the t — J model away from half filling is 
possible at present, we focus on the so called spin-fermion model. This Hamiltonian can 
be derived in fourth order degenerate perturbation theoryiE from the p — d, three band, 
Emery modelil, that gives a detailed description of the cuprate materials. There the role 
of perturbation is played by the hybridization term between the p-orbital (oxygen) and the 
d-orbital (copper). A further simplification of the model was proposed by Zhang and Ricell, 
that leads to the t — J model. 

A typical copper-oxide two leg ladder, as those present in Sr 14 _ a .Ca a; Cu2404i is depicted 
in Fig. [I]. It is generally accepted that the dopant holes reside on p-orbitals on the oxygens 
sites, whereas on the Cu 2+ ions a localized hole resides, represented by spin 1/2 operator 
which interact via a nearest neighbor exchange. 

The spin-fermion Hamiltonian is defined as follows: 



H SF = t < ')" ' • -^X>-S, 

+J H Y,Si-Si>- (1) 

(*,*') 

The index i (j) runs over the Cu (O) sites, cj CT creates a hole in an oxygen p band and 
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FIG. 1: Schematic picture of a two leg ladder Copper-Oxide. 

Sj are spin operators for the copper ions. The coefficients aijj take care of the sign of the 
p-d overlap and a^j = 1 if j = i + \x or i + \y and otij = 2 if j = i — |x or i — \y. Finally 
the operator Rj is defined as 



Ri = E (-^'^ia^C^. (2) 

{j,k;i),i,a,/3 

Followingil we can define the following operator centered on the copper site Pj )(T = 
(1/2) Yl(j-i) (~ l) Qi ' J Cj i(T which represents non orthogonal orbitals with a high weight on the 
i site. Their anti-commutation relations are 

and we can rewrite the Hamiltonian in terms of these operators as follows 



1=1. ..L 1=1.. .L 

A=1,2,ct A=l,2,a,/3 

,A • S J+1| a, (4) 

1=1. ..L 1=1. ..L 

A=l,2 



L is the number of the rungs along the ladder and A = 1, 2 distinguishes the two legs. For 
the sake of generality, an anisotropy in the Heisenberg term is allowed. 

The different steps of our procedure are the following: first find orthogonal (Wannier 
states) for the holes, then go to a (coherent states) path integral formulation for spins and 
fermions and perform the Gaussian integration of the fermionic degrees of freedom. The 
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remaining part of the calculation is devoted to the evaluation of the resulting Fermi deter- 
minant in the long- wavelength low-frequency limit. This expansion includes the coupling 
constant Jk to all order. 

Wannier states are easily find via Pt )CT = a/ e(k)/k jCr where e (k) — (l — COh ( fc3;a )+^ os fe a )/ 2 ^ _ 
Here a is the lattice constant and we used a two dimensional Fourier transform where k y 
takes only values and it /a distinguishing between symmetric (bonding) and antisymmetric 
(antibonding) states. The partition function can be expressed as a path integral 

Z = J D[f*]D[f]D[n]e- s ^ : (5) 
where Ssf = Sh + S s . The action S s contains all terms with spins degree of freedoms onlyEI: 



/ft ✓ ) ( ) 

dr -iS A (6i,a) • + Hu cis (SQ (t) 

J, A 



(6) 



where Cl is a unimodular field, S is the spin per site (1/2 in our case) and A is the vector 
potential for a (Dirac) monopole: e abc (dA a / 'd&j,) = Cl c . 

It is by now well accepted that the effective low energy field theory of the <i-dimensional 
Heisenberg antiferromagnetic model is given by the (d + 1) NLaMEHl. In the case of a 
ladder one obtains the (1 + 1) NL<tM00. For this reason, here we will deal mainly with the 
part of the action which contains fermionic degrees of freedom Sh'- 



S h=^2 ^*.« ( iuJn + Ate ^ ~ V) 5 k, q 5 a ,f3 + gy/e(k)e(q)(T Qi p ■ Ci k _ q f q ^ } (7) 

kqa.fi 

here k = (k x , k y ,u n ) where u n = n (2n + 1) / (3 are the fermionic Matsubara frequency and 
g = 4J K S. It is natural to decompose the inverse propagator into G^ 1 = G$ 1 — £ where 
the free part is 



G = (iu n + 4te(k) - /i) 5 Kq 5, 



a,/3, 



(8) 



and the fluctuating external potential is 



s = -5 , V e ( k )e(q)cr a , / 3 • VL k _ q . 



(9) 
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Since, according to Eq. (|7|) the action Ssf is bilinear in the fermionic variables, we can 
integrate them out. This leads to Ssf = S s — trlnG -1 . Defining the matrix 

A = ^/7^)8 k , q 8 a , p (10) 

and a rescaled propagator G~ l through 

G- 1 = A^G-'A- 1 (11) 

we can write 

tr In (G- 1 ) = tr In (AA f ) + tr In (g~*) , (12) 

the first term gives just a constant and we can ignore it. Again we decompose the rescaled 
inverse propagator as G^ 1 = Gq 1 — X which brings us to 

A-i f iu> n + Ate Ck) — _ _-, , , , , 

G Q = I --^j-j j OkqOa/3 = g {K x ,k, y ,U n )0k,q0 a fi, (13) 

E = -gtlk-^Lo-v ■ cr a/3- (14) 

The remaining part of the calculation is devoted to the evaluation of Sh e fr = — tr In (^G^ 1 ^ 
in the continuum limit. 

A. Parameterizations 

As we already mentioned, in the undoped regime where no holes are present, it has 
proven very effective a mapping from a antiferromagnetic Heisenberg spin ladder to a (1+1) 
NLaM. This mapping rely on the idea that although long range order (here antiferromag- 
netic) is prohibited in one dimension, the most important contribution to the action are 
given by paths in which antiferromagnetic order survives at short distance. Accordingly the 
dynamical unimodular field is decomposed in a Neel modulated field n plus a ferromagnetic 
fluctuating contribution. A gradient expansion in the dynamical field brings then to the 
(1+1) NLaM. The gradient expansion is justified when the correlation length of the spin 
is much larger than the lattice constant a. However the prediction of the NLcrM, i. e. a 
finite correlation length and a triplet of massive modes above the ground state0'@i!0 re- 
main valid until £ ~ 2.5a as numerical calculations on the isotropic Heisenberg ladder have 
shownlll. 
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The basic assumption of this work is then that such a parameterization is still meaningful 
as long as the spin liquid state is not destroyed by doping, as seems to be the case in 



n iy \ and 1^ are two slowly varying, orthogonal, vector fields describing locally anti- 
ferromagnetic and ferromagnetic configurations, respectively, is normalized such that 
l n i,A| 2 = 1- The lattice constant a in front of l^x in eq. ([15|) makes explicit the fact that 1^ 
is proportional to a generator of rotations of n^, namely to a first-order derivative of n^. 

In the particular geometry of a ladder, this decomposition give rise to two local order 
parameters, n^i and n^- However we assume that spins across the chain are rather strongly 
correlated such that they will sum up to give rise to an antiferromagnetic configuration, or 
subtract and give a ferromagnetic fluctuation. A further parameterization is then 



with Nj ■ Mi = and |N;| 2 = 1. 

The next step is the gradient expansion, or equivalently, in Fourier space, an expansion 
in powers of k. In (1+1) dimensions the field N will get no scaling dimension, whereas 
the fields 1 and M get scaling dimension -1. Accordingly, in the subsequent expansion we 
will need to keep terms with up to two derivative and any power of the field N. Terms 
containing 1, M are marginal whenever two fields or one field and one derivative are present. 
Higher order terms are irrelevant and will be discarded. This correspond to expand all our 
quantities up to O (a 2 ). 

The self energy has then the following expansion 



experiments, where a finite spin-gap is also seen in the doped casd'Hl'Sll. Then, 
ref.0, we parameterize the spin field in the following way 



as e.g. m 




(15) 




(16) 



E = E 00 + Eoi + E 02 + Ei + E 2 + O (a 3 ) , 



(17) 
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where the various quantity are 



500 — —g8k y -q y ,'ir'Nk x -q a .+n,ui- v -(T cl( p, (18) 

5 1 = —ag5k y -q y ,o'M-k x -q x +TT,u,-v-0'a{3, (19) 

s o 2 = ^K-^ (N M 2 ) kx _ qx+7r ^ ■ a aP , (20) 

51 = ^-lk-q,£j-y - &aP, (21) 

5 2 = ^(N|l| 2 ) k _ q+Qi _-0" a/3 , (22) 

where Q = (ir/a,7r/a) is the antiferromagnetic modulation vector suitable for a ladder 
geometry. We also regroup the zero-th order term in F ^ = Gq — Sqo- 

The evaluation of the various contribution in the continuum limit, proceeds very similarly 
as in ref.El, and we refer to that paper for a more detailed explanation. The quantity to be 
evaluated is 

S h eff = -tr In (F- 1 ) - trln(l — F (S i + S 02 + Si + S 2 )) . (23) 
We need then to find the inverse of F^ 1 up to O (a). It turns out that 

F = FD- 1 - aFD^RD- 1 + O (a 2 ) , (24) 
where the various matrices are 

D = D(k,uj)8 kq 8 aP , (26) 

R = -g5 ky - qyt7T ^2 ( k r - Qr + 8 r , x 7r/a) drg^ik^^k-g+Q ■ a a/3 , (27) 

r=x,T 

and we used the shorthand notation 

g x (k, u n ) = g x (k + Q, cu n ) , (28) 

D(k,u n ) = g l (k, uj n ) g l (k, uj n ) - g 2 . (29) 

We first consider the term 

tr In (F' 1 ) = tr In (G x ) + tr In (l - G Soo) • (30) 
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The second term of this equation is reduced to the calculation of 

oo _ 

£-tr((Wf, (31) 

m=l 

where each term has the following expansion 

tr (Go^oo) = (g) m ^2 9a(k)go(k + ofigoik + q-z)ga(k + q 4 ) ■ ■ ■ g (k + qm-i)g (k + q m ) 

k,q2...q m 

xn*;n;% 3 • • • N^tr •••<), (32) 

with m an even integer. The trace over the Pauli matrices can be carried out using a trace 
reduction formulaS The gradient expansion in Eq. (|32D is then obtained by performing 
an expansion of the product of propagators go(k) ■ ■ -go(k + q m ) in powers of the variables 
q 2 , (fe, . . . q m that appear as argument of the vector field N. The result obtained is@ 



trln(l-G Soo 
with the definition 



dxdr 



^ |^ N | 2 + |d r N| 2 



(33) 



Xa/3 



d 2 



dq a dq p 



In [l-g 2 g(k)g(k + q + Q)] 5 qy0 



(34) 

This does not present 



We can now pass to the evaluation of the second term in Eq. 
particular problems, since after expanding all the quantities, it reduces to the evaluation of 
a finite number of traces. The result is 

trln (1 — F (£01 + S 02 + £1 + £ 2 )) = Jdxdr X r (N x «9 T N) ■ (l x + 1 2 ) 

<r f dxdrx(h + h) 2 ■ (35) 
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Here we omitted to write a Gaussian term oc M 2 , completely decoupled, which can be 
integrated out without further consequences. The quantities \t and x are given by 



X 



-i D-Wa^toWD-^k + Q), 

k 

(g 1 (k + Q)-g \k))} 2 . 



(36) 
(37) 



They are generalized susceptibilities of the holes in presence of long-wavelength spin 
fields. In particular the zeros of D(k) determine the dispersion of such holes. The bands 
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FIG. 2: Effective hoies lowest-band emerging from our theory. Parameters are t = 0.24, Jk = 1 eV.The 
minimum falls exactly at (ak) = 2tt/3 

originating in such a way correspond to free holes moving in a staggered magnetic field. 
Such a staggered field would break translation invariance by one site and we would obtain 
four bands in the reduced Brillouin zone. Instead in our procedure we never broke explicitly 
translation invariance, so that we obtain genuinely two bands in the Brillouin zone. The 
lowest of these two band is symmetric in character (bonding). In Fig. ||] we show it for values 
of the constants relevant for the Copper-Oxide ladder i. e. a band-width of ~ 0.5 evil and 
This band is in good agreement with accurate calculations on the one hole 
spectrum of the t — J model. In particular, in the isotropic t — J model, for t/J~2 the 
same qualitative feature are observed: a global maximum at (ka) = 0, global minima at 
(ka) « ±27r/3 and local maxima at (ka) = ±7r0Il. 

Now that we calculated the long wavelength contribution coming from the holes, we still 
have to consider the continuum limit (in the low energy sector) of the pure spin action S s 
given by eq. (|G|). The result is 



The very last step is the Gaussian integration of the l± field, leaving us with the effective 
long-wavelength action for the antiferromagnetic order parameter, a (1+1) NLcrM: 




dx dr (li + 1 2 ) 2 



(38) 
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Sett — Sh efr + S s eg — — / dxdr 



v\d x N\ 2 + -\d T N\ 2 



(39) 



where the NLcrM parameters are given by 

/ 



./ 



g2j Xxx 



V 



4J||+2J ± + §^ 



IT 



(40) 



v = a 



Xxx \ 

2 I 



[4J||+2Ji 



(41) 



Hence, the spin-fermion model with mobile holes interacting with an antiferromagnetic 
background is mapped into an effective NLcrM whose coupling constant depend on doping 
through the generalized susceptibilities in Eqs. (|34f) , fl36D , and Q37|) . 

Now we can immediately transpose to our model of a doped spin liquid, some known 
result for the NLcrM, e. g. mainly the presence of a gap which separates the singlet ground 
state from a triplet of magnetic excitations. This gap should persist as long as the continuum 
approximation is valid. 

The fact that the NLcrM in (1+1) dimension has a gap above the ground state can be 
established in a variety of ways. Using the two loop beta functioni§ one obtains 



A 



_2w I 2~K 



(42) 



where A is a cutoff of the order of the inverse lattice constant. Now we have an explicit 
analytic form for the doping dependence of the spin gap in the spin-liquid state of a two leg 
ladder. 

To study the behavior of the gap with doping we have to distinguish two regimes where 
the lowest effective band has minimum either at zero or at 2/3n. For Jk > 2t the minima fall 
in ±2/37r. Here all the generalized susceptibilities in Eqs. (@), (|36|) , and ( j37D contribute to 
lower / and, since from eq. fl42|) A is an increasing function of /, they make the gap smaller 
for any value of the constants (see Fig. ^). This is comforting, since, as we mentioned, for 
Jk very large the physics of the Spin-Fermion model should be similar to that of the t — J 
modelil, and for that one, TDMRG simulations show that the gap decreases at least in a 
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FIG. 3: J K > 2t. (a) Generalized susceptibilities of Eqs. (|34|), <^§), and (|7p for J K =2,t = 0.76 eV. For 
Jk > 2t all the susceptibilities contribute to lower / hence the gap decreases for small doping for any value 



of the constants, (b) Normalized gap of eq. (42). Here we fixed the exchange constants to J|| = Jj_ = 0.108 
eV. A is the gap without doping. 




an 



FIG. 4: J K < 2t. (a) Generalized susceptibilities of Eqs. (§J), ©, and @ for J K = 3t = 1.8 eV. For 
Jk < 2i, one susceptibility, \, grows with doping and contributes to increase / and hence the gap. For 
(Jk, t) ^S> ( J\\ , J±) we can have an increasing gap for small doping, (b) Normalized gap of eq. (|42|). Fixing 
the exchange constants to J|| = Jj_ = 0.108 eV is enough to have an increasing gap for small doping. 

strong anisotropic case (Jj_ = 10J||). When Jx < 2t the band minimum falls in zero and 
there is one susceptibility, x, which instead makes / grow. /. In this regime there is then a 
(small) region of parameters where the gap grows with doping (see figure f§). 

Before passing to a comparison with experiments, we want to comment on a possible 
simplifying understanding. 

A simple picture to explain the observed diminishing of the spin gap with doping in 
Sr^-zCa^C^/iO/u, is that (at least for low doping concentration where speaking of a spin 
liquid is still feasible) the effect of the holes is that of renormalizing the anisotropy parameter 
A = J±/J\\ for the spin part towards larger values. In many studies on the 2 leg ladder 
Heisenberg antiferromagnet000, the spin gap is seen to increase with A. In fact, the same 
occurs in the NLcxM without doping in the range A ~ 1 -r- 2. 

We can now pass to our mapping of a doping spin liquid to an effective NLcxM. According 
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to equations ( f40| , f4TD effective coupling constants J\\ , J± can be denned for the doped system 
such that the form of the NLaM parameters is that for a pure spin system@ i. e. 



f = l shk> (43) 



2oSJ||,/l + i. (44) 
A small doping expansion in the regime Jk > 2t leads to 

J„ = J, + \^j^l 5 + 0(5*), (45) 



J i = J i - 



(!(£^ + 2(4J|| + 2Jl) + MLt^), + 0( , 2) , (46) 



so indeed Jy, Jj_ are seen respectively to increase, decrease, such that A decreases. However, 
such an interpretation breaks down beyond 5 ~ 0.04 whereas /, v are still well defined 
positive constants. This means that beyond such doping, this simplified picture cannot be 
naively applied and holes have a more effective way of lowering the gap. 



III. COMPARISON WITH EXPERIMENTS 



We come now to the comparison with experiments. Our theory depends on four pa- 
rameters t, Jk, J\\, J± which we now want to fix to physical values. ARPES experiment on 
Sri4Cu2404i were performed by Takahashi et. al. H who found a band matching the period- 
icity of the ladder with a bandwidth of ~ 0.5 -f- 0.4 eV. Adjusting our lowest band to have 
such a bandwidth we obtain a relation between t and Jk- On the other hand, experiments 
on the Cu02 cell materials and band theory calculation!!©!! agreed in assuming a value of 
Jk of the order of Jk ~ 1 -r- 2 eV. This in turn gives us a value of t ~ 0.24 0.76 eV, which 
is also consistent with the same calculation. 

The debate around an anisotropy of the spin exchange constantsHCI in Sri4_ x Ca :r Cu2404i 
seems now to be resolved in favor of isotropy or light anisotropy of the coupling constant: 
J±/J\\ ~ 0.8@. We adjusted the value of the momentum cutoff A by fixing the theoretical 
gap with the experimental one for the undoped compound Sri 4 Cu2404i. Finally, to compare 
with the measured values of the gap for different doping concentration x in Sr 14 _ a .A x Cu 2 404i 
(where A can be either divalent Ca 2+ ,Ba 2+ or trivalent Y 3+ ,La 3+ ), we still need a relation 
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between the A substitution x and the number of holes per copper site present in the ladder 
5. This is another unsettled issue of the telephone number compound. In particular Osa- 
fune et. al. I studying the optical conductivity spectrum, inferred that with increasing Ca 
substitution x, holes are transferred from the chain to the ladder. On the other hand Niicker 
et. al. @ argue that in the series compound Sr 14 _ 2 .Ca x Cu2404i the number of holes in the 
ladder is almost insensitive to Ca substitution x (although a small increase is observed). 
Here we will assume that Sr^^A-zC^/tO/u is an example of doped spin liquid and will use 
the data from!. The result of our theory can be seen in figure [|. There we used isotropic 
exchange constant, but the theoretical curve did not change in a visible way if an anisotropy 
of J±/ J\\ ~ 0.8 was inserted. We see from the figure that the spin gap becomes zero for 
5 ~ 0.37, beyond this value the coupling constants / and v would become imaginary signal- 
ing that our effective model cease to make sense. This means that for such doping ratios 
our parameterization fljjD is no longer valid, in the sense that it does not incorporate the 
most important spin configurations. However our theory could cease to make sense much 
before. If one takes the point of view of the t — J model (as we said the Spin-Fermion model 
should map to it for large Jk) the holes introduced in the system couple rigidly to the spins 
forming singlet with the P$ states. In the worst case this would limit the correlation length 
of the spin to the mean hole-hole distance 1/5. In our case this happens at a doping ratio 
of 5 » 0.15. 

A word of caution should be mentioned with respect to comparison with experimental 
results. A still unresolved controversy is present between NMFiiS! and neutron scatter- 
ingii experiments, where the latter see essentially no doping dependence of the spin gap. 
Without being able to resolve this issue, we would like, however, to stress, that beyond the 
uncertainties in experiments, the doping behavior obtained for the spin-gap agrees with the 
numerical results in TDMRG and is opposite to the one obtained in mean-field treatments, 
making clear the relevance of fluctuations. 

IV. CONCLUSION 

In this paper we studied the behavior of the spin gap of a two leg Heisenberg antifer- 
romagnetic ladder as microscopically many holes are introduced in the system. Such a 
situation can be physically realized in the series compound Sri4_. E A. E Cu2404i with A=Ca, 
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1000 , , , , , , , , , , , , , , , r 




holes per copper 5 



FIG. 5: Result of our theory and comparison with experiments. The values of the constants used in equation 
( fj2|) are t = 0.76, Jk = 2 eV, J\\/J± = 1. The momentum cutoff A was fixed by fixing the the value of the 
gap with the one measured in Sr^Cu^O^. For anisotropic case J_\_/J\\ = 0.8 the curve does not change 
appreciably. 

Y, La, and numerous result are now available from experiments. On the theoretical side, 
however, there is a contradiction between previous analytical treatments on the one hand, 
and TDMRG simulations or NMR experiments on the other hand. Whereas in the first 
case, a magnon gap increasing with doping is predicted, a decrease is observed in accurate 
numerical simulation and experiments. 

Starting from the spin-fermion model we were able to solve the contradiction using a con- 
trolled analytical treatment that properly takes into account fluctuations in the continuum 
limit. Integrating out the fermions we were left with a Fermi-determinant which we can 
evaluate exactly in that limit. The result is a non linear a model with doping dependent 
parameters. The spontaneously generated mass gap of this theory is seen to decrease as 
holes are introduced. Once physical value for the parameters are given, we obtained very 
good agreement with NMR experiments performed on Sri4_ x A x Cu2404i. 
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